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Abstract B037 is of interest because it is both the most luminous and the most highly red- 
dened cluster known in M3 1 . Images of deep observations and of highly spatial resolutions 
with the Advanced Camera for Surveys on the Hubble Space Telescope (HST) firstly show 
that this cluster is crossed by a dust lane. Photometric data in the F606W and F814W filters 
obtained in this paper provide that, colors of (F606W - F814W) in the dust lane are red- 
der ~ 0.4 mags than ones in the other regions of B037. The HST images show that, this 
dust lane seems to be contained in B037, not from the M3 1 disk or the Milky Way. As we 
know, the formation of dust requires gas with a rather high metallicity. However, B037 has 
a low metallicity to be [Fe/H] - -1.07 + 0.20. So, it seems improbable that the observed 
dust lane is physically associated with B037. It is clear that the origin of this dust lane is 
worthy of future study. In addition, based on these images, we present the precise variation 
of ellipticity and position angle, and of surface brightness profile, and determine the struc- 
tural parameters of B037 by fitting a single-mass isotropic King model. In the F606W filter, 
we derive the best-fitting scale radius, ro - 0.56 + 0.02" (- 2.16 ± 0.08 pc), a tidal radius, 
r, = 8.6 + 0.4"(= 33. 1 + 1. 5 pc), and a concentration index c = log(r,/ro) = 1.19+0.02. In the 
F814W filter, we derive ro = 0.56 + 0.01" 2.16 + 0.04pc), r, = 8.9 ±0.3" 34.3 ±1. 2 pc), 
and c - log(r,/r()) = 1.20 + 0.01. The extinction-corrected central surface brightness is yuo - 
13.53+0.03 mag arcsec^^ intheF606W filter, and 12.85+0.03 mag arcsec^^ in the F814W fil- 
ter, respectively. We also calculate the half-light radius, at r/, - 1.05 ±0.03"(= 4.04 + 0.12pc) 
in the F606W filter and r/, = 1 .07 + 0.01"(= 4.12 + 0.04 pc) in the F8 14W filter, respectively. 
In addition, we derived the whole magnitudes of B037 in V and / bands by transforming the 
magnitudes from the ACS system to the standard system, which are in very agreement with 
the previous ground-based broad-band photometry. 

Key words: galaxies: evolution - galaxies: individual (M31) - globular cluster: individual 
(B037) 



1 INTRODUCTION 

Globular clusters (GCs) are effective laboratories for studying stellar evolution and stellar dynamics, and 
they are ancient building blocks of galaxies which can help us to understand the formation and evolution 
of their parent galaxies. In addition, GCs exhibit surprisingly uniform properties, suggesting a common 
formation mechanism. 

The closest other populous GC system beyond the halo of our Galaxy is that of M3 1. The study of 
M3 1 has been and continues to be a keys tone of extraga lactic astronomy jBarmbv et al.l l2000l) . and the 
study of GCs in M31 can be traced back to|Hubbld (ll932|). M3 1 GC B327 rB for 'Baade ' ) or B o37 (Bo for 
'Bologna', see Battistini 1987), which, in the nomenclature introduced by Huchra et alj ( 199 lb is referred 
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to as B037, a designation from the Revised Bologna Catalogue (RBC) of M31 GCs and candidates (Galleti 
et al. 2004, 2006, 2007), which is the main catalo g used in studies of M31 GCs. The extremely red color of 
B037 was firstl y noted by Kron & Mavalll ((1960), who suggested that this cluster must be highly reddened. 
Two years later. lVetesnikId 1962a) de termined magnitu des of 257 M31 GC candidates including B037 in the 
UBV photometric system, and then lVetesnikI (Il962bl) studied the intrinsic colors of M31 GCs, and found 
that B037 was the most highly reddened with E(B-V) - 1 .28 in his sample of M3 1 GC candid ates based on 
the photometric catalog of Vetesnik ( 1962a). With low -resolution spectroscopv. lCrampton et al. (1985) also 
found that B037 is the most highly reddened GC candida te in M3 1 to have ^(B-y) = 1.48. Based on a large 
database of multicolor photometry. lBarmby et alJ (l2000 i) determined the reddening value for each individual 
M31 GC including B037 using the correlations between optical and infrared colors and metallicity by 
defining various "reddening-free" parameters, an d the reddening valu e of B037 is E(B - V) - 1.38 + 0.02 
(which is kindly given us by P. Barmby). Again, Barmbv et alj (l2002b ) derived the reddening value for this 
cluster to be E(B - V) - 1.30 + 0.04, using the spectroscopic metallicity to predict the intrinsic colors. 
iMa et alj (j2006a) also determine the reddening of B037 by comparing independently obtained multicolor 
photometry with theoretical stellar population synthesis models to be E(B - V) = 1.360 + 0.013, which is 
in good agreement with the other results. Following the methods of iBarmbv et al.l ( l2000l) . lFan et alj ( |2008l) 
(re-)determined reddening va lues for 443 clusters and cluster candidates including B037, and the redding 
value of B037 obtained bv .Fan et all (I200&) is E{B - V) = 1.21 + 0.03, which is a little smaller than the 
previous determinations. 

The brightest GCs i n M3 1 are more luminous than the most brightest Galactic cluster, oj Centauri. 
Among these are B037 jvan den Berghi [1968) and Gl (see details from Barmbv et al.. ,2002b) . These 
two clusters are both considered as the possible remna nt core of a former dwa r f galaxy whi c h lost 
most of its enve lope through tidal interactions with M31 jMevlan & Heggiel Il997t iMevlan et all 1200 ll: 
iMackev & vand en Bergh, 2005; Ma et al., 2006b, 200l). 

In this paper, we will present the photometric data of B037 using its deep images obtained with the 
Advanced Camera for Survey (ACS) on the HST. The deep images of highly spatial resolutions showed 
that this cluster is crossed by a dust lane. Our results provide that colors of F814W-F606W in the dust lane 
are redder ~ 0.4 mags than ones of the other regions. In addition, we studied structures of B037 in detail 
based on these images. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 HST images of B037 

We searched the HST archive and found B037 to have been observed with the ACS-Wide Field Channel 
(WFC) in the F606W and the F814W filters, which were observed on 2004 August 2 and on 2004 July 4, 
respectively. The exposure time is 2370.0 seconds for both bands. The HST ACS-WFC resolution is 0.05" 
per pixel. The images in F606W and F814W both show that B037 is crossed by a dust lane. Fig. 1 clearly 
shows the dust lane, which crosses B037. If the dust lane is true, its color should be different from ones of 

the other regions. 

If not otherwise stated, the magnitudes are always on the VEGAMAG scale as defined by lSirianni et al.l 

(l2005h . The relevant zero-point for this system is 26.398 and 25.501 for WFC F606W and WFC F814W, 
respectively. A distance to M31 of 780 kpc (1" subtends 3.85 pc) is adopted in this paper 

2.2 Color difference between the dust lane and the other regions 

In order to study whether the color difference between the dust lane and the other regions in B037 exists, 
we select nine points, three of which (No. 7, 8 and 9) are located in the dust lane, the other six are randomly 
located in the other regions (see Figure 2). For each sample point, the PHOT routine in DAOPHOT ( Stetson, 
Il987h is used to obtain magnitude. We adopt an aperture of a diameter of 4 pixels. The photometric data 
for these nine sample points are given in Table 1, in conjunction with the Icr magnitude uncertainties 
from DAOPHOT. Column 4 gives the color of (F606W - F814W). From Table 1, we can see that colors of 
(F606W - F814W) in the dust lane are redder ~ 0.4 mags than ones of the other regions. 
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Fig. 1 The images of GC B037 observed in the F606W and F814W filters of ACS/HST. The 
images clearly show that the cluster is crossed by a dust lane. The image size is 17.5" x 17.5" 
for each panel. 
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Fig. 2 The sample positions of photometry (black circles) are showed in the image of GC B037 
observed in the F606W filter of ACS/HST. An aperture of radii of 2 pixels is adopted for pho- 
tometry. The image size is 17.5" x 17.5". 



Table 1 Photometric data for B037 



Source 


F606W 


F814W 


F606W-F814W 


No. 


(mag) 


(mag) 


(mag) 


1 


23.29 + 0.26 


21.35 ±0.16 


1.94 


2 


23.03 + 0.23 


21.14 ±0.15 


1.89 


3 


22.95 + 0.22 


21.00 ±0.14 


1.95 


4 


23.22 + 0.25 


21.25 ±0.15 


1.97 


5 


23.17+0.25 


21.36 ±0.16 


1.81 


6 


21.15 ±0.10 


19.10 ±0.06 


2.05 


7 


22.93 ± 0.22 


20.53 ±0.11 


2.40 


8 


23.41 ±0.28 


20.94 ±0.13 


2.47 


9 


24.66 ± 0.49 


22.31 ±0.25 


2.35 



Detailed study of B327 based on HST images 



5 



2.3 Surface brightness profiles 

We used the iraf task ellipse to obtain F606W and F8 14W surface brightness profiles for B037. B037 center 
position was fixed at a value derived by object locator of ellipse task, however an initial center position was 
determined by centroiding. Elliptical isophotes were fitted to the data, with no sigma clipping. We ran two 
passes of ellipse task, the first pass was run in the usual way, with ellipticity and position angle allowed 
to vary with the isophote semimajor axis. In the second pass, surface brightness profiles on fixed, zero- 
ellipticity isophotes were measured, sin ce we choose to fit circ ular models for the intrinsic cluster structure 
and the point spread function (PSF) as iBarmbv et al. I ( l2007h did (see §2.4 for details). The background 
value was derived as the mean of a region of 100 x 100 pixels in "empty" areas far away from the cluster. 

2.3.1 Ellipticity and position angle 

Tables 2 and 3 give the ellipticity, e = 1 - b/a, and the position angle (P.A.) as a function of the semi- 
major axis length, a, from the center of annulus in the F606W and F814W filter bands, respectively. These 
observables have also been plotted in Figures 3 and 4, respectively; the errors were generated by the iraf 
task ELLIPSE, in which the ellipticity errors are obtained from the internal errors in the harmonic fit, after 
removal of the first and second fitted harmonics. From Table 3, and Figs. 3 and 4, we can see that, the 
values of ellipticity and positio n angle ca nnot be obtained within 0.1448" in the F814W filter because 
of very high elliptic ity (> 1.0). iMa et alj (l2006b.) analyzed the same F606W image of B037 used here, 
fitting a lKind (1 19621) model to a surface brightness profile made from a PSF-deconvolved image. They also 
plotted the distributions of ellipticity and the position angle as a function of the semi-major axis length. 
Comparison of Fig. 2 of Ma et al. (2006b) and Figs. 3 and 4 sho ws that, the genera l trend of the cluster's 
ellipticity as a function of semimajor axis radius is similar between lMa et alj ( l2006bl) and the present paper. 
The comparison also shows that uncertainties in the exact value of the PA are only of secondary importance 
for the general trend in ellipticity observed, given that the PA determination between Ma et al. (2006b) and 
the present paper differs somewhat greatly. There are a numb er of possible r easons for the offsets in PA 
observed between these two studies. The main reason is that. iMa et all (l2006b ) used the PSF-deconvolved 
image. Other reasons include those related to the positions of the centering of isophotes and the different 
geometrical parameters set when fitting. In addition. Fig. 3 shows that the ellipticity varies significantly 
with position along the semimajor axis radius, especially smaller than 0.5". In the F814W filter band, the 
ellipticity is larger than 1.0 along the semimajor axis radius smaller than 0.1448". 

2.4 Point spread function 

At a distance of 780 kpc, the ACS/WFC has a scale of 0.05 arcsec - 0.19 pc pixel"', and thus M31 clusters 
are clearly resolved with it. Their observed core structures, however, are still affected by the PSF. We chose 
not to deconvolve the d ata, instead fitti ng structural models after convolving them with a simple analytic 
description of the PSF as lBarmbv et"an (,2007.) did. To estimate the PSF for the WFC. lBarmbv et al. I ( l2007h 
used the iraf task ellipse with circular symmetry enforced to produce intensity profiles out to radii of about 
2" (40 pixels) for a number of isolated stars on a number of images, and combined them to produce a single, 
average PSF. This was done separately for the F606W and F814W filters. They originally tried to fit these 
with simple Moffat profiles (with backgrounds added), but found that a better description was given by a 
function of the form below. For the combination of the WFC and F606W filter, 

/psF = /o [l + (-R/a.'0686)^]"' '' , (1) 

which has a full width at half-maximum of FWHM = 0'.'125, or about 2.5 px; for the combination of the 
WFC and F814W filter, 

/psF = /o [l + (-R/a.'0783)^]"' '' , (2) 

which has a full width at half-maximum of FWHM - 0'.'14 5, or about 2.9 px. In addition, since this PSF 
formula is radially symmetric and the models of King (il966h we fit are intrinsically spherical, the convolved 
models to be fitted to the data are also circularly symmetric. 
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Table 2 B037: Ellipticity, e, and position angle (P.A.) as a function of the semimajor axis, a, in 
the F606W filter of HST ACS-WFC 



a 

(arcsec) 




6 


P.A. 

(deg) 


a 

(arcsec) 




6 


PA. 

(deg) 


0.0260 


0.638 


± 0.228 


92.9 ± 15.5 


0.3757 


0.177 


±0.031 


69.8 ± 


5.6 


0.0287 


0.638 


± 0.229 


93.2 ± 15.6 


0.4132 


0.151 


± 0.029 


64.7 ± 


5.9 


0.0315 


0.639 


± 0.230 


93.4 ± 15.7 


0.4545 


0.090 


± 0.027 


60.3 ±9.1 


0.0347 


0.640 


± 0.232 


93.7 ± 15.8 


0.5000 


0.005 


± 0.025 


172.6 ± 


30.0 


0.0381 


0.642 


± 0.233 


94.0 ± 15.9 


0.5500 


0.060 


± 0.020 


155.5 ± 


9.8 


0.0420 


0.643 


± 0.235 


94.4 ± 16.0 


0.6050 


0.117 


±0.015 


156.3 ± 


3.9 


0.0461 


0.645 


± 0.236 


94.7 ± 16.0 


0.6655 


0.174 


±0.012 


157.2 ± 


2.2 


0.0508 


0.647 


+ 0.182 


95.2 ± 12.3 


0.7321 


0.233 


±0.011 


157.2 ± 


1.5 


0.0558 


0.599 


±0.159 


96.8 ± 11.2 


0.8053 


0.278 


±0.011 


159.1 ± 


1.3 


0.0614 


0.546 


±0.142 


98.5 ± 10.6 


0.8858 


0.322 


±0.011 


160.8 ± 


1.1 


0.0676 


0.503 


±0.127 


100.2 ± 10.0 


0.9744 


0.358 


±0.012 


162.1 ± 


1.2 


0.0743 


0.458 


± 0.099 


102.2 ± 8.3 


1.0718 


0.380 


± 0.021 


164.5 ± 


2.1 


0.0818 


0.400 


± 0.059 


104.3 ± 5.5 


1.1790 


0.367 


± 0.022 


168.0 ± 


2.2 


0.0899 


0.410 


± 0.050 


101.0 ±4.6 


1.2969 


0.343 


± 0.025 


169.2 ± 


2.7 


n nr» on 


n /no 
U.420 


± U.U44 


no "7 1 Q n 
98. / ± 3.9 


1 .4ZDO 


n T 1 n 


± U.UZJ 


166.8 ± 


2.8 


0.1088 


0.437 


± 0.046 


97.6 ± 3.9 


1.5692 


0.252 


± 0.022 


165.1 ± 


3.0 


0.1197 


0.428 


± 0.028 


96.6 ± 2.5 


1.7261 


0.239 


± 0.021 


165.7 ± 


2.9 


0.1317 


0.410 


± 0.027 


96.5 ± 2.5 


1.8987 


0.211 


± 0.026 


163.6 ± 


4.0 


0.1448 


0.400 


±0.031 


96.3 ± 3.0 


2.0886 


0.201 


± 0.029 


152.8 ± 


4.6 


0.1593 


0.364 


± 0.023 


95.1 ±2.3 


2.2975 


0.188 


± 0.037 


150.1 ± 


6.3 


0.1752 


0.352 


± 0.027 


94.8 ± 2.8 


2.5272 


0.182 


± 0.033 


149.9 ± 


5.8 


0.1928 


0.337 


± 0.027 


93.2 ± 2.9 


2.7800 


0.180 


± 0.034 


145.6 ± 


6.0 


0.2120 


0.311 


± 0.027 


92.5 ± 3.0 


3.0580 


0.191 


±0.031 


137.5 ± 


5.2 


0.2333 


0.287 


± 0.026 


90.6 ±3.1 


3.3638 


0.143 


± 0.034 


125.4 ± 


7.2 


0.2566 


0.258 


± 0.026 


88.6 ±3.3 


3.7001 


0.180 


± 0.041 


121.1 ± 


7.1 


0.2822 


0.233 


± 0.027 


85.4 ±3.8 


4.0701 


0.257 


± 0.033 


121.6 ± 


4.1 


0.3105 


0.207 


± 0.029 


81.2 ±4.5 


4.4772 


0.233 


± 0.048 


121.6 ± 


6.5 


0.3415 


0.189 


± 0.030 


75.3 ±5.0 


4.9249 


0.237 


± 0.063 


116.1 ± 


8.5 



2.5 Extinction 

When we fit models to the brightness p rofiles of B037, we will correct the inferred magnitude parameters 
for extinction. The reddening law from Cardelli et al.l (Il989 ) is employed in th is paper. The effectiv e wave- 
lengths of the ACS F606W and F814W filters are A^g = 5918 and 8060 A (Sirianni et al.', '2005*), so that 
from Cardelli et al. (1989), Af606w - 2.8 x£'(B - V) and Apguw - 1.8x£'(B-y)(seeBarmbv et al. ,2007; 
iMcLaughUn et all I2008L for details). The reddening value of E(B - V) = 1.360 + 0.013 from Ma et al] 
(l2006ah is adopted in this paper 



2.6 Magnitudes of B037 in F606W and F814W filters 

We derived the total flux of B037 in F606W and F814W filter bands using the iraf task phot in dapphot as 
below: measuring aperture magnitudes in concentric apertures with an interval of 0.1", drawing magnitude 
growth curves, and paying attention to where the flux does not increase. At last, we obtained the magnitudes 
of B037 in F606W and F814W to be 16.21 + 0.010 and 14.16 + 0.006, respectively. In the photometry, we 
derived the background value as the mean of a region far away from the cluster (see §2.3 for details). 
We use VEGAMAG photometric system. In order to allow a meaningful comparison with the previous 
ground-based broad-band photometry of Barmby et al. (2000), we transformed the magnitudes from the 
ACS system to the standard broad-band photome tric system by following the transformation equations and 
coefficients o f Table 22 of Isirianni et aL 1 (l2005h . The results are mv(ACS) - 16.83 (this paper) versus 
my = 16.82 dBarmbv et al.ll2000l) , and m^A CS) = 14.15 ( t his pa per) versus ;«/ = 14.16 dBarmbv et all 
l2000h . Our results are in good agreement with lBarmbv et alj (l2000b . 



Detailed study of B327 based on HST images 



7 



Table 3 B037: Ellipticity, e, and position angle (P.A.) as a function of the semimajor axis, a, in 
the F8 14W filter of HST ACS-WFC 



a 

(arcsec) 




6 


P.A. 

(deg) 


a 

(arcsec) 




e 


PA. 

(deg) 


0.0260 








0.4132 


0.031 


± 0.030 


77.8 ± 28.3 


0.0287 








0.4545 


0.031 


± 0.029 


20.2 ± 27.5 


0.0315 








0.5000 


0.044 


± 0.027 


161.2 ± 17.6 


0.0347 








0.5500 


0.084 


± 0.023 


150.3 ± 8.3 


0.0381 








0.6050 


0.127 


±0.021 


150.3 ±5.0 


0.0420 








0.6655 


0.175 


±0.019 


153.6 ±3.4 


0.0461 








0.7321 


0.220 


±0.016 


157.4 ± 2.3 


0.0508 








0.8053 


0.247 


±0.013 


161.8 ± 1.7 


0.0558 








0.8858 


0.251 


±0.014 


167.8 ± 1.8 


0.0614 








0.9744 


0.263 


±0.017 


170.0 ±2.1 


0.0676 








1.0718 


0.293 


± 0.034 


170.8 ±4.0 


0.0743 








1.1790 


0.297 


± 0.035 


172.5 ±4.1 


0.0818 








1.2969 


0.230 


± 0.028 


171.5 ±4.0 


0.0899 








1.4266 


0.216 


± 0.025 


166.7 ± 3.7 


0.0989 








1.5692 


0.198 


±0.031 


165.5 ±5.1 


0.1088 








1.7261 


0.198 


± 0.025 


169.5 ± 4.0 


0.1197 








1.8987 


0.188 


± 0.029 


167.5 ±5.0 


0.1317 








2.0886 


0.139 


±0.031 


166.9 ± 6.9 


0.1448 








2.2975 


0.117 


±0.031 


117.8 ± 8.0 


0.1593 


0.908 


±0.117 


89.9 ± 6.9 


2.5272 


0.100 


± 0.034 


148.1 ± 10.2 


0.1752 


0.878 


± 0.026 


90.9 ± 1.5 


2.7800 


0.118 


±0.051 


141.0 ± 13.3 


0.1928 


0.827 


±0.151 


90.8 ± 9.6 


3.0580 


0.094 


± 0.043 


115.8 ± 13.7 


0.2120 


0.749 


± 0.034 


90.3 ± 2.3 


3.3638 


0.094 


± 0.035 


132.7 ± 11.1 


0.2333 


0.731 


± 0.037 


89.4 ± 2.6 


3.7001 


0.103 


± 0.056 


121.3 ± 16.4 


0.2566 


0.695 


±0.041 


87.1 ±2.9 


4.0701 


0.127 


±0.061 


120.1 ± 14.6 


0.2822 


0.624 


±0.031 


84.4 ± 2.3 


4.4772 


0.162 


±0.031 


115.7 ±5.9 


0.3105 


0.546 


± 0.035 


80.4 ± 2.6 


4.9249 


0.091 


± 0.054 


131.4 ± 17.7 


0.3415 


0.401 


± 0.035 


72.4 ± 3.2 


5.4174 


0.150 


± 0.065 


135.5 ± 13.3 


0.3757 


0.258 


± 0.029 


63.4 ±3.8 


5.9591 


0.188 


± 0.044 


161.8 ±7.4 



3 MODELS AND FITS 
3.1 Structural models 

After elliptical galaxies, GCs are the best understood and most thoroughly modelled class of stellar systems. 
For example, a large majority of the ~ 150 Galactic GCs have been fitte d b y the s i mple models of single- 
mass, isotropic, lowered isothermal spheres developed bv lMichiS (Il963h and lKing|(ll966l) (hereafter "King 
models"), yielding comprehensive catalogs of cluster structural parameters and physical properties (see 
iMcLaughlin & van der Marel, 2005, and references therein). For extr agalactic GCs, HST imaging data have 
been used to fit King models to a large number of GCs in M 3 1 (e.g.. lBarmbv et al.L l2002al: iBarmbv et al. 



2007', and references the rein), in M33 dLarsen et all l2002b . and in NGC 5128 (e.g.. [Harris etall 12002 : 
McLaughl in et al.l |2008[ and references therein). In this paper, we fit the usual King models to the density 
profile of B037 observed with ACSAVFC. 



3.2 Observed data 

Tables 4 and 5 list the surface brightness, fi, of B037, and its integrated magnitude, m, as a function of radius 
in the F606W and F814W filters, respectively. The errors in the surface brightness were also generated by 
the iRAF task ellipse, in which they are obtained directly from the root mean square scatter of the intensity 
data along the zero-ellipticity isophotes. In addition, the surface photometries at radii where the ellipticity 
and position angle cannot be measured, are obtained based on the last ellipticity and position angle as the 
IRAF task ELLIPSE is designed. 
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Fig. 3 Ellipticity as a function of the semimajor axis in the F606W and F814W filters of 
ACS/HST. 



3.3 Fits 

Our fitting procedure involves computing in full large numbers of King structural models, spanning a w ide 
range of fixed values of the appropriate shape parameter Wo (see iMcLaughlin & van der Mare]Ll2005L m 
detail). And then the models are convolved with the ACS/WFC PSF for the F606W and F814W filters of 
equations of (1) and (2); 



Imod(R'/ro) X /psF [(x - x'), (y - y')] dx' dy' , 



(3) 



where /^od = /mod/^o; and /psp is the PSF profile normalized to unit total luminosity (see lMcLaughlin et all 
I2OO8, in detail). We changed the luminosity density to surface brightness Ji*^^^ - -2.5 log W^^^] before 
fitting them to the observed surface-brightness profile of B037, fj. - - 2.5 log [/(/?/ro)//o], finding the 
radial scale ro and central surface brightness yUo which minimize for every given value of Wq. The 
(Wo, ro,fio) combination that yields the global minimum over the grid used defines the best-fit model 
of that type: 



X 



(4) 
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in which o", is the error in the surface brightness. Estimates of the one-sigma uncertainties on these basic fit 
parameters follow from their extreme values over the subgrid of fits withx^/v < X^iJ'^ + 1' here v is the 
number of free parameters. Figure 5 shows our best King fits to B037. In Fig. 5, open squares are ellipse 
data points included in the least-squares model fitting, and the asterisks are points not used to constrain 
the fit. These observed data points shown by asterisks are inc luded in the radius of R < 2 pixels = (X'l, 
and the isophotal intensity is dependent on its neighbors. As Barmbv e t al. I (12007 ') pointed out that, the 
ELLIPSE output contains brightnesses for 15 radii inside 2 pixel, but they are all measured from the same 13 
central pixels and are not statistically independent. So, to avoid excessive weighting of the central regions 
of B037 in the fits, we only used intensities at radii Rmin, Rmin + (0.5, 1 .0, 2.0 pixels), or R > 2.5 pixels as 
iBarmby eT al. (2007) used. Table 6 summarizes the results obtained in this paper. 



3.4 Comparison to previous results 



iMa et all (l2006b ) analyzed the same F606W image of B037 used here, fitting a iKingl (1 19621) model to a 
surface brightness profile der ived from a PSF-d econvolved image. They derived the scale radius ro - 0'.'72 
(it is called the core radius in lMa et all (l2006bl) ). half-light radius ri, = I'.'l 1, concentration index c - 0.91, 
and central surface brightness yu(0) - 17.21 mag arcsec"^ (using the v alue for ex tinction adopted in this 
paper, this becomes /do = 13.40 mag arcsec"^). Comparing the results of iMa et al] (l2006b;) with Table 6 of 
this paper, we find that our model fits produce a somewh at higher con cen tration and s maller scale radius. 
These differences come from: (i) using different models dKind (1 19621) vs lKingI (Il966l) ). (ii) the observed 
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Fig. 4 P. A. as a function of the semimajor axis in the F606W and F814W filters of ACS/HST. 
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Fig. 5 Surface brightness profile of B037 measured in the F606W and F814 filters. Dashed 
curve (blue) trace the PSF intensity profiles and solid (red) curves are the PSF-convolved best-fit 
models. Open squares are ellips data points included in the;^'^ model fitting, and the asterisks are 
points not used to constrain the fits (see the text in detail). 



Detailed study of B327 based on HST images 



11 



Table 4 B037: Surface brightness, ji, and integrated magnitude, m, as a function of the radius in 
the F606W filter of HST ACS-WFC 



R 








m 


R 










m 


(arcsec) 




(mag) 


(mag) 


(arcsec) 




(mag) 




(mag) 


0.0260 


17 


327 ± 


0.007 


23.827 


0.3757 


17 


792 ± 





040 


18.456 


0.0287 


17 


328 ± 


0.008 


23.827 


0.4132 


17 


863 ± 





046 


18.264 


0.0315 


17 


328 ± 


0.008 


23.827 


0.4545 


17 


944 ± 





049 


18.123 


0.0347 


17 


329 + 


0.009 


23.827 


0.5000 


18 


040 ± 





049 


17.962 


0.0381 


17 


330 ± 


0.010 


23.827 


0.5500 


18 


148 ± 





048 


17.827 


0.0420 


17 


331 ± 


0.011 


23.827 


0.6050 


18 


267 ± 





048 


17.672 


0.0461 


17 


331 ± 


0.012 


23.827 


0.6655 


18 


389 ± 





053 


17.549 


0.0508 


17 


332 ± 


0.014 


22.086 


0.7321 


18 


495 ± 





061 


17.414 


0.0558 


17 


334 ± 


0.015 


22.086 


0.8053 


18 


598 ± 





070 


17.295 


0.0614 


17 


338 ± 


0.016 


22.086 


0.8858 


18 


716 ± 





077 


17.165 


0.0676 


17 


341 ± 


0.018 


22.086 


0.9744 


18 


854 + 





079 


17.039 


0.0743 


17 


346 ± 


0.020 


21.452 


1.0718 


19 


006 ± 





077 


16.927 


0.0818 


17 


351 ± 


0.022 


21.452 


1.1790 


19 


193 ± 





107 


16.822 


0.0899 


17 


356 ± 


0.024 


21.452 


1.2969 


19 


440 ± 





105 


16.731 


0.09S9 


1 / 


363 ± 


0.027 


IlAjZ 


1.4/66 


19 


721 ± 





1 16 


16.642 


0.1088 


17 


372 ± 


0.027 


21.062 


1.5692 


20 


001 ± 





113 


16.570 


0.1197 


17 


382 ± 


0.029 


20.552 


1.7261 


20 


293 ± 





116 


16.505 


0.1317 


17 


395 ± 


0.028 


20.552 


1.8987 


20 


597 ± 





122 


16.451 


0.1448 


17 


411 ± 


0.028 


20.370 


2.0886 


20 


872 ± 





128 


16.401 


0.1593 


17 


429 ± 


0.029 


19.964 


2.2975 


21 


224 ± 





105 


16.358 


0.1752 


17 


450 ± 


0.028 


19.964 


2.5272 


21 


457 ± 





112 


16.320 


0.1928 


17 


475 ± 


0.026 


19.764 


2.7800 


21 


719± 





138 


16.283 


0.2120 


17 


504 ± 


0.026 


19.528 


3.0580 


22 


082 ± 





154 


16.251 


0.2333 


17 


538 ± 


0.025 


19.337 


3.3638 


22 


603 ± 





164 


16.225 


0.2566 


17 


578 ± 


0.026 


19.091 


3.7001 


23 


042 ± 





225 


16.206 


0.2822 


17 


624 ± 


0.026 


19.009 


4.0701 


23 


694 ± 





467 


16.191 


0.3105 


17 


675 ± 


0.030 


18.802 


4.4772 


24 


509 ± 





571 


16.182 


0.3415 


17 


732 ± 


0.034 


18.634 


4.9249 


25 


173 ± 


1 


342 


16.172 



data are obtained with different ways . In ("iil.lMa et al ] (l2006bh derived the surface brightness profile from a 
PSF-deconvolved image; in addition. 'M a et al.l (12006b') derived the surface brightness profile with ellipticity 
and position angle allowed to vary with the isophote semimajor axis, however, in this paper, we derived the 
surface brightness profile on fixed, zero- ellipticity isophotes, s ince we choose to fit circular models for the 
intrinsic cluster structure and the PSF as Barmbv et al. ] (l2007b did (see §2.4 for details). In fact, from Fig. 
5 of this paper and Fig. 3 of Ma et al. (2006b), we can see that, the observed data are somewhat different 
between Ma et al. |(j2006b) and this paper. 

iBarmbv et al. I f2007) analyzed the same F606W and F8 14W images of B037 used here with nearly the 
same o bserved data and method. The results of comparison are listed in Table 7 (Table 5 of iBarmbv et al. I 
(120071) in the electronic edition did not list the results of B037 in the F814 W filter.), from which we can see 
that the results obtained in this paper are in good agreement with ones of IBarmbv et al. I (l2007b (about the 
central surface brightnesses, we have corrected them using the value for extinction adopted in this paper). 



4 DISCUSSION AND SUMMARY 



As discussed in §3.1, it is impossible that the dust lane comes from the MiUcy Way. Another possibility 
is that the dust lane is contained in B037 itself. As we know, the formation of dust requires gas with a 
rather high metallicitv. lPerrett et al.l (l2002h presented metallicities for more than 200 GCs in M3 1 including 
B037, using the Wide Field Fibre Optic Spectrograph at the 4.2 m William Herschel Telescope in La Palma, 
Canary Islands, which provides a total spectral coverage of ~ 3700-5600 A with two gratings. One grating 
(H2400B 2400 line) yielded a dispersion of 0.8 A pixel"' and a spectral resolution of 2.5 A over the range 
3700-4500 A covering the CN feature at 3883 A, the H and K lines of calcium, H6, the CH G band and the 
4000 A continuum break, and the other grating (R1200R 1200 line) presented a dispersion of 1 .5 A pixel ' 
and a spectral resolution of 5.1 A over the range 4400-5600 A to add absorption features such as HyS, 
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Table 5 B037: Surface brightness, ju, and integrated magnitude, m, as a function of the radius in 
the F8 14W fiher of HST ACS-WFC 



R 








m 


R 








m 


(arcsec) 




(mag) 


(mag) 


(arcsec) 




(mag) 


(mag) 


0.0260 


15 


301 ± 


0.010 


21.800 


04132 


15 


772 ± 


0.032 


16.190 


0.0287 


15 


302 ± 


0.011 


21.800 


0.4545 


15 


863 ± 


0.036 


16.048 


0.0315 


15 


303 ± 


0.012 


21.800 


0.5000 


15 


967 ± 


0.038 


15.889 


0.0347 


15 


303 + 


0.013 


21.800 


0.5500 


16 


078 ± 


0.037 


15.754 


0.0381 


15 


304 ± 


0.015 


21.800 


0.6050 


16 


190 ± 


0.033 


15.598 


0.0420 


15 


305 ± 


0.016 


21.800 


0.6655 


16 


300 ± 


0.036 


15.474 


0.0461 


15 


306 ± 


0.018 


21.800 


0.7321 


16 


407 ± 


0.046 


15.338 


0.0508 


15 


308 ± 


0.019 


20.061 


0.8053 


16 


543 ± 


0.053 


15.218 


0.0558 


15 


310 ± 


0.021 


20.061 


0.8858 


16 


706 ± 


0.054 


15.094 


0.0614 


15 


313 ± 


0.024 


20.061 


0.9744 


16 


847 ± 


0.054 


14.974 


0.0676 


15 


317 ± 


0.026 


20.061 


1.0718 


16 


995 + 


0.061 


14.868 


0.0743 


15 


321 ± 


0.030 


19.430 


1.1790 


17 


185 ± 


0.105 


14.767 


0.0818 


15 


326 ± 


0.033 


19.430 


1.2969 


17 


459 ± 


0.086 


14.681 


0.0899 


15 


331 ± 


0.037 


19.430 


1.4266 


17 


717 ± 


0.107 


14.596 


0.0989 


15 


337 ± 


0.041 


19.430 


1.5692 


18 


006 ± 


0.101 


14.527 


n 1 noo 
U. lUSo 


1 ^ 


347 ± 


0.043 


ly.Uiy 


1. //Ol 


18 


279 ± 


0.102 


14.464 


0.1197 


15 


360 ± 


0.043 


18.531 


1.8987 


18 


567 ± 


0.086 


14.411 


0.1317 


15 


367 ± 


0.045 


18.531 


2.0886 


18 


833 ± 


0.092 


14.359 


0.1448 


15 


375 ± 


0.048 


18.350 


2.2975 


19 


167 ± 


0.095 


14.316 


0.1593 


15 


396 ± 


0.046 


17.940 


2.5272 


19 


460 ± 


0.094 


14.277 


0.1752 


15 


411 ± 


0.046 


17.940 


2.7800 


19 


649 ± 


0.156 


14.240 


0.1928 


15 


428 ± 


0.045 


17.738 


3.0580 


20 


075 ± 


0.137 


14.207 


0.2120 


15 


447 ± 


0.046 


17.496 


3.3638 


20 


538 ± 


0.103 


14.181 


0.2333 


15 


466 ± 


0.048 


17.301 


3.7001 


21 


002 ± 


0.138 


14.160 


0.2566 


15 


501 ± 


0.046 


17.046 


4.0701 


21 


399 ± 


0.203 


14.142 


0.2822 


15 


531 ± 


0.043 


16.960 


4.4772 


21 


964 ± 


0.228 


14.129 


0.3105 


15 


580 ± 


0.037 


16.745 


4.9249 


22 


519± 


0.240 


14.117 


0.3415 


15 


632 ± 


0.031 


16.573 


54174 


23 


311 ± 


0.588 


14.106 


0.3757 


15 


695 ± 


0.029 


16.388 


5.9591 


23 


508 ± 


0.782 


14.096 



Table 6 Structural parameters of B037 



Parameters F606W F814W 

To 0.56 ±0.02" (= 2.16 ± 0.08 pc) 0.56 + 0.01" (= 2.16 ± 0.04 pc) 
r, 8.6 + 0.4" (= 33.1 ± 1.5 pc) 8.9 + 0.3" (= 34.3 + 1.2 pc) 
c = log(r,/ro) 1.19 + 0.02 1.20 ±0.01 
n, 1.05 ±0.03"(= 4.04 ± 0.12 pc) 1.07 ± 0.01"(= 4.12 ± 0.04 pc) 
fip (mag ai-csec-^) 13.53 ± 0.03 12.85 ±0.03 



Table 7 Results of comparison 



F606W 


F814W 


Parameters 


Barmbvetal. (2007) 


This paper 


Barmbvetal. (2007) 


This paper 




ff.'56 


0'.'56 


a' 59 


a'56 


c = log(r,/ro) 


1.23 


1.19 


1.18 


1.20 




1'.'09 


I'.'OS 




1'.'07 


yUo (mag arcsec"^) 


13.45 


13.53 


12.75 


12.85 



the Mg b triplet, and two iron line s near 5300 A. Then.lPerrett et alJ (l2002h calculated 12 absorption-line 
indices based on the prescription of lBrodie & Huchral(ll99ol) . By the comparison of the Une indices \yith the 
published M3 1 GC [Fe/H] values from the previous literature (iBonoh et al.lll987llBrodie & Huchrall 19901: 
iBarmbv et al.l 12000) . the results of linear fits were obtained. Final cluster metalUcities were determined 
from an unweighted mean of the [Fe/H] values calculated from the CH (G), Mg b, and Fe53 Une strengths. 
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For B037.| Perrett et alJ (l2002h obtained its metallicity to be [Fe/H] = - 1 .07 + 0.20. It is clear that B037 has 
a low metallicity. So, it is intricate that where is the dust lane from? 

In this paper, using the images of deep observations and of highly spatial resolutions with the ACS/ HST, 
we firstly present that the GC B037 in M3 1 is crossed by a dust lane. Photometric data in the F606W and 
F814W bands provide that, colors of (F606W - F814W) in the dust lane are redder ~ 0.4 mags than ones 
in the other regions of B037. From the HST images, this dust lane seems to be contained in B037, not from 
the MiUcy Way. However, the formation of dust requires gas with a rather high metallicity. So, it seems 
impossible that the observed dust l ane is physically ass ociated with B037 itself, which has a low metallicity 
to be [Fe/H] = -1 .07 + 0.20 froml Perrett et all (l2002h . So, that the observed dust lane in the view of B037 
is from B037 itself or from the Milky Way needs to be confirmed in the future. In addition, based on these 
images, we present the precise variation of ellipticity and position angle, and of surface brightness profile, 
and determine the structural parameters of B037 by fitting a single-mass isotropic King model. 
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